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We study epidemic processes with immunization on very large 1-dimensional lattices, where at 
least some of the infections are non-local, with rates decaying as power laws p(x) ~ x~ tT ~ 1 for 
large distances x. When starting with a single infected site, the cluster of infected sites stays 
always bounded if a > 1 (and dies with probability 1, if its size is allowed to fluctuate down to 
zero), but the process can lead to an infinite epidemic for a < 1. For a < the behavior is 
essentially of mean field type, but for < a < 1 the behavior is non-trivial, both for the critical 
and for supercritical cases. For critical epidemics we confirm a previous prediction that the critical 
exponents controlling the correlation time and the correlation length are simply related to each 
other, and we verify detailed field theoretic predictions for a \ 1/3. For a = 1 we find generic 
power laws with continuously varying exponents even in the supercritical case, and confirm in detail 
the predicted Kosterlitz-Thouless nature of the transition. Finally, the mass N(t) of supercritical 
clusters grows for < a < 1 like a stretched exponential. This implies that networks embedded 
in 1-d space with power-behaved link distributions have infinite intrinsic dimension (based on the 
graph distance), but are not small world. 



I. INTRODUCTION 

Epidemic spreading, both on regular lattices and on 
random graphs, have attracted increasing attention in 
the statistical physics community [TJ [2]- Basically one 
distinguishes between two types of epidemics. In both 
all agents except for the seeds are susceptible, while in- 
fectious agents (such as the seeds in the initial configu- 
ration) stay infectious for a finite time (during one time 
step in the following), which means that they can infect 
agents that share with them a link. After this infective 
period, they either become susceptible again (SIS model), 
or they become "removed" , i.e. either immune or dead 
(SIR model). In both cases, agents stay at fixed places 
(whence they can be identified with lattice sites or net- 
work vertices), so that the entire dynamics is contained 
in the changes of their internal (S, I, R) states. 

There exist of course several generalizations of these 
two basic schemes (e.g. cooperative infection [3JE], in- 
fection using the "power of choice" 5 , moving agents \§\ , 
memory effects [7], etc.), but even the basic schemes show 
rich behavior, if we allow for different network topologies. 

In the following we shall only deal with basic SIR epi- 
demics, with discrete time and infective period equal to 
one time step. In this case, the process must die out on 
any finite system, since susceptibles are used up but not 
replenished, and the set of "removed" sites becomes, for 
large times, just a percolation cluster. More precisely, we 
shall discuss two models, both of which involve infection 
over large distance, with the infection probability decay- 
ing as some inverse power of the distance 8 10 . In both, 
the sites can be viewed as sites in an infinitely large 1-d 
lattice (in the simulations we use L = 2 64 sites, which is 
big enough so that we never have to worry about finite 
size effects, except in a few cases pointed out later |11|). 
In both models, each infected site first attempts to in- 
fect ko other sites (and succeeds so with probability 1, if 
these sites are still susceptible), and then attempts with 



probability p to infect one more site. Thus, the average 
number of newly attempted infections per node, which 
is also equal to the average out-degree of the generated 
graph, is 



fcout = ko + P ■ 



(1) 



But both models differ slightly in how these sites are 
connected, i.e. how the attempted infections are chosen: 

• In model (A) we assume that in each attempt the 
target site is chosen randomly, with a distance ±x 
from the infectious site that is distributed according 
to a power law, 



P(x) 



(2) 



for large x. More precisely, this distribution is ob- 
tained by first drawing a real-valued random num- 
ber y uniformly from the interval (1/L CT ,1] (for 
a < 0), [1, L a ) (for a < 0) or [0,lnL] (for u = 0) 
and setting 



. _ f L?r 1/CT J for a + 0, 
" S \e y \ for a = 0. 



(3) 



Notice that we do not check that all attempts try 
to infect different targets. If several attempts are 
aimed at the same target, all except the first one 
are simply lost. 

• In model (B) we first infect the left and right neigh- 
bors, and only in the subsequent infections distant 
sites are chosen, again with the same probability 
P(x) given above. 

In model (A) , it can happen that all attempts during 
one time step try to infect target sites that are no longer 
susceptible, in which case the epidemic dies. This cannot 
happen in model (B). There, the right neighbor of the 



2 



rightmost infected site is always susceptible, as is also the 
left neighbor of the leftmost site. Model (B) is indeed, 
as far as the geometric structures obtained for t — » oo 
are concerned, a modification of the Watts-Strogatz [TJ] 
small world model that was studied previously in [T5HT%] . 

In the next section, we will discuss the most important 
features common to both models. In Sec. Ill we shall 
treat in more detail the case a = 1, where we find a 
number of non-trivial exact results. The critical case of 
model (A) with < a < 1 is studied in Sec. IV [model 
(B) is always supercritical, as it never can die, thus no 
critical phase exists] . The supercritical case for < a < 1 
is finally discussed in Sec. V. The paper concludes with 
a discussion in Sec. VI. 



II. GENERAL FEATURES 



s? 10 



10 



10 



15 



10"' 



10° 



10 

x 



10 



1 


1 1 

t=1 - 




t=2 - 




t-3 




t=4 




t=5 - 




t=6 




t 7 











10 



15 



We always start with a single seed located at the origin, 
x = 0. Boundary conditions are periodic. The jump 
probability P(x) is cut off at x = L, i.e. the length of 
the widest jumps allowed is precisely the lattice size. The 
distribution of infected (or "active" , as we shall call them 
in the following) sites at time t > is denoted as p(x,t). 
The distribution of removed sites is then 

t-i 

R(x,t) = J2p( x ,t')- ( 4 ) 

t'=0 

The average number of active sites at time t is then 

n(i)=J>(s,t), (5) 

X 

while the number of immune sites is 

t-i 

AT(t) = £ii(M) = £n(t'). (6) 

x t'=0 

Notice that N(t)+n(t) can also be interpreted as the av- 
erage number of sites reached by chains of at most t links 
from a randomly chosen pivot. Therefore, if N(t) ~ t D , 
the exponent D would be the intrinsic (or "topological" 
19J [20 ) dimension of the final network. 

For model (B), R(x,t) = 1 and p{x,t) = for all 
|x| < t. Thus the active sites form two outgoing waves, 
one moving to the right and the other to the left. For 
a > the infection is sufficiently short ranged that these 
two waves don't interfere with each other, in the limit 
t oo. Basically the same is also true for model (A), 
although in that case R(x,t) and p(x,t) are no longer 
strictly 1 and 0, respectively, for \x\ < t. Nevertheless, 
activity dies out also there for any finite x, i.e. p(x, t) — > 
for t — > oo. 

This scenario no longer holds for a < (for the case 
a = see [l5j). In that case, P{x) would not be nor- 
malizable for infinitely large lattices, hence finite lattice 
effects become important. Indeed, on large but finite lat- 
tices, both R(x,t) and p{x,t) become uniform, see Fig.JT] 



FIG. 1. (Color online) Densities of active sites at times t = 
1, 2, ... 7 for model (A) with a = —0.1, ko = 2, and p = 1, i.e. 
with exactly 3 offsprings per active site. Because of left-right 
symmetry, we show here and in the following plots only the 
distribution for x > 0. The densities become more and more 
uniform for increasing t, showing the mean field nature when 
a < 0. 
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FIG. 2. (Color online) Densities of active sites at times t = 
5, 10, . . . 75 for model (A) with a = 0.5, ko = 1, and p = 1, 
i.e. with exactly 2 offsprings per active site. Notice the tails 
decaying x'"" 1 for all times. For small x the activity dies out 
because (i) the density of immune sites gets higher and higher, 
and (ii) the parent activity shifts more and more outside. 

Also, the formation of finite loops is suppressed by pow- 
ers of L (the number of loops with m links, each limited 
in length to less than some constant, scales as L 1+ma ), 
so that the cluster of immunes becomes locally tree like 
in the limit L — > oo and the model becomes of mean field 
type. 

In contrast, for a > it is found that both R(x,t) 
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FIG. 3. (Color online) Densities of immune sites at at the 
same times as in Fig. [2] 



and p(x, t) decay asymptotically as \x\ a 1 for all times 
t > 0, see e.g. Figs. [2| [| (see also [2T]). This is easily 
understood. First of all, they cannot decay faster, since 
the offspring distribution of a population concentrated at 
the origin would decay like that and any smearing due to 
a finite extend of the parent population can only make the 
offspring distribution wider. On the other hand, if p(x, t) 
decays at some given time t not slower than |x| _<7_1 , then 
the distribution of its offsprings is given for sufficiently 
large \x\ by 



(7) 



In this expression we have neglected saturation effects 
(not all infections are successful, because not all sites are 
susceptible) , but this approximation should be correct for 
large \x\ where most of the sites are susceptible. Since 
R(x,t) is just a sum over p(x,t') with t' < t, it satisfies 
the same asymptotic behavior. 

Notice that this argument only tells about the limit 
where we first let x — ¥ oo, and then let t become large. 
It does not prove that R oc (x) — lim^oo R(x, t) and 
Poo(x) = Ihnt-j.oo p( x jt) decay asymptotically as |a;|~ cr ~ 1 . 
For < a < 1 this seems to be correct nevertheless (see 
Figs. 2j3), but it it docs not hold for a > 1. In that 
case the process dies with probability one for model (A) 
(for an elegant and simple proof, see [22] ) , but it survives 
forever for model (B). In the latter case the average num- 
ber of active sites, n{t) = p(x, t), tends to a constant 
and the wave of active sites has a stationary profile in a 
co-moving frame (in a frame moving with constant veloc- 
ity the profile widens due to fluctuations of the velocity) . 
Profiles R(x,t) are shown in Fig. [4] for model (B) with 
a = 1.5 and one long-range contact per site. We see that 



FIG. 4. (Color online) Densities of the immune cluster at 
times t = 10, 20, 40, . . . 1280 for model (B) with a = 1.5, ko = 
3, and p = 1, i.e. with two local and one long-range offsprings. 
For large x, all densities decay ~ l/x rT+1 , but for large times 
this tail sets in later and later. 
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this behavior sets in later and later for increasing t. The 
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FIG. 5. (Color online) Plot of log log n(t) versus \ogt, where 
n(t) is the average number of active sites, for the same run 
shown also in Figs. [2] and [3] The straight line represents 
a stretched exponential fit to the data with 50 < t < 75, 
corresponding to 2.6 x 10 6 < n(t) < 1.6 x 10 8 . Statistical 
errors are much smaller than the symbols. 



bulk of the outgoing wave has finite width (i.e., becomes 
increasingly sharper in a log- log plot such as Fig. [4] ). 
In the terminology of [53], the fronts for model (B) are 
pushed when a > 1. 

The same is true for model (A), although one has to 
condition on (exponentially rare) surviving events in or- 
der to see this. This could be done by using e.g. PERM 
[2"1] . but we have not done it. 

Notice that the times in Fig. [4] are exponentially in- 
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creasing. The linear progression of the front in the log- 
log plot then means that the immune cluster travels at 
constant speed. This is in contrast to the case a < 1 
shown in Figs. [2] and [3] (where time increases linearly 
between successive curves). There the wave of infection 
travels with a speed that increases faster than a power 
with time, in agreement with [211 125] but in contrast to 
[21)]). On the other hand, the decreasing distances be- 
tween successive curves in Figs. [2] and [3] show that this 
increase of speed is less than exponential, in contrast to 
[2"Tl |2"5] . In order to see this more clearly, we plot n(t) 
in Fig. [5] for the runs shown in Figs. [2] [3] More pre- 
cisely, Fig. 5] shows loglogn(i) plotted against logt. An 
exponential n{t) ~ exp(at) would correspond in this plot 
to a straight line with slope 1. This is obviously not 
observed (statistical errors are smaller than the symbol 
sizes). Rather, the data for large t suggest a stretched 
exponential n(t) ~ exp(at 7 ) with 7 = 0.59(1). But we 
should be careful in accepting this fit as the true asymp- 
totic scaling. First of all, the data in Fig. [2] are slightly 
curved, even for the largest i, suggesting that this esti- 
mate of 7 is too high. Also, fitting stretched exponentials 
is notoriously fraught with uncertainties. The same be- 
havior is also seen for model (B) (data will be shown 
later). It shows also that the mass of clusters with di- 
ameter £ increases slower than exponentially with £, in 
contrast to claims made in |16| (exponential increase is 
of course seen in the mean field regime, a < 0). More 
details will be given in Sec. V. 



III. THE CASE cr = 1 

It is well known that the case of interacting Levy flights 
with a = 1 is very special, in particular in one dimension 
of space. This was first found by Dobrushin, Ruelle [27] . 
and Dyson [28] who showed that 1-d Ising models with 
long range interactions can only have a finite temperature 
phase transition, if a < 1. Very soon after this, Anderson 
et al. [29] and Thouless [30] showed that for cr = 1 one not 
only does have a phase transition, but that this transition 
is similar to the Berezinskii-Kosterlitz-Thouless (BKT) 
transition in displaying a region with generic power laws 
in the supercritical phase with cr-dependent exponents. 
As in the XY model, the reason is that configurations 
can be described as organized by defects which inter- 
act with each other by an attractive logarithmic poten- 
tial. This argument was later extended by Cardy [31] to 
Potts and other models. Due to the Fortuin-Kasteleyn 
[52"] relationship between the Potts model and percola- 
tion [51] , this applies also to percolation and thus also to 
SIR epidemics. But it seems that the consequences for 
the latter have never been worked out in detail, with one 
notable exception: It was shown in 33J that the perco- 
lation transition for cr = 1 is discontinuous in the sense 
that the order parameter (the density of infected sites 
for t — > 00) jumps discontinuously, when p is increased 
through the percolation threshold. This might seem con- 
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FIG. 6. (Color online) Log-log plot of n(t) versus t for model 
(A) with a = 1, and for different values of fc ou t- All curves for 
fcout > 2.6 seem to become asymptotically straight, giving a 
rough estimate of fc out = 2.6(1) for the critical point. 
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FIG. 7. (Color online) Data collapse for model (A), us 
ing the finite-t scaling ansatz Eq. (111. Parameters 
k c = 2.602, r\ = 1.42, t = 0.88, 
points with t > 3 are plotted. 



and y — 1.37. 



are 
Only 



tradictory to the claim of universal power laws (which 
usually hold only at continuous phase transitions), but 
several similar "hybrid" cases, where aspects typical of 
a first order transition coexist with aspects of a second 
order transition, have been found recently also in other 
contexts [31153]. 

Numerical results for the increase of n(t) in both mod- 
els at cr = 1 are shown in Figs. [6] to [8] From Fig. [6] we 
see that model (A) exhibits indeed generic power law be- 
havior (in agreement with the predictions of [2UM31] ) for 
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FIG. 8. (Color online) Exponents rj(k out ,a = 1) governing 
the increase of active sites obtained from Figs. [6] and [9] For 
both models, 77 is plotted against the number of long-range 
outgoing links, i.e. it is plotted against k out for model (A) 
and against k out — 2 for model (B). 



all fc ou t > k c , where k c w 2.6(1): 



(8) 



At the critical point, a straightforward fit gives i] = 
r?(fc out ) = 1.42(4). 

More precise determinations of k c and rj are possible 
by using the finite (cluster-)size scaling expected for BKT 
transitions. Let us define e — fc ou t — k c . Near the critical 
point there exists a characteristic time scale r(e) which 
diverges as e — > 0. Inversely, we can for each t define an 
effective distance from the critical point as e(t). Using the 
latter we can make the finite-time scaling (FTS) ansatz 



n(t,k out ) ~t"g[e/e(t))], 



(9) 



where g(z) is an analytic function joining smoothly the 
sub- and supercritical regions. While this ansatz is com- 
mon to models with short and long range infections, the 
scaling of e(t) in the limit t — > oo is different. For short 
range contacts it is a power law, while for a BKT transi- 
tion we expect 



e(t) ~ [In 



to 



Equations ^ and ( 10 ) can be combined to 



n(t, k a 



t v g[{k 



*c)(m r n 

to 



(10) 



(11) 



A data collapse based on Eq. pdj ) is shown in Fig. [7j We 
see a nearly perfect collapse (only points with t > 3 are 
plotted), giving our best estimates 



k c = 2.60(1), 77 = 1.42(1), t = 0.9(1) and y 



1.37(3). 

(12) 



TABLE I. Exponents ^(fc ou t) for model (B). Column 2 gives 
our results, column 3 is from Table 1 of [18], noticing that 
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1 + 77 and 13 = (fc out - 2)/2. 
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To our knowledge, neither 77 nor y have been calculated 
before. We conjecture that they are universal for all mod- 
els where P(x) ~ 1/x 2 asymptotically and where, in con- 
trast to model (B), the epidemic can die (we made also 
preliminary simulations of a generalization of model (B) 
where left and right neighbors are infected with probabil- 
ities 0.9. The results support the conjecture). Our value 
of k c is consistent with the exact bound k c > 2 for 
this class of models. 

Combining Eqs. Q and (11) gives that g(z) behaves 
for large z as a stretched exponential, In g(z) <~ z 1 ^ , and 
that 



r](kout) - V ~ (fcout - k c ) 1/y 



(13) 



This is reasonably well satisfied. 

As seen from Fig. |9j essentially the same behavior is 
found also for model (B), with one important exception: 
Since model (B) with a > 1 is supercritical for all fc ou t > 
2, all curves become straight lines for t — > 00 and r](k out ) 
tends to zero for fc out — > 2. The dependence of ri(k ou t) 
on fc ou t is shown in Fig. [8] and Table 1. The fact that 
n(t, fcout) increases as a power of t is indeed known [TH 
118) . In |15j exact upper bounds on 77(fc ou t) were given for 
small fc out , and these were compared to simulation results 
in [18] . When comparing our estimates with the results 
of |18j . we should notice that I + 77 is the graph dimension 
of the cluster of immunes, and that the constant f3 used 
in [TH] corresponds to (fc ou t — 2)/2. From Table 1 we 
see that our data are roughly 10 times more precise, but 
otherwise they are in very good agreement. 

For both models 77 becomes the same for large fc ou t) 
which is to be expected: For large fc ou t, the evolution 
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FIG. 9. (Color online) Same as Fig. [6] but for model (B). 
Since model (B) is supercritical for all fc out > 2, all curves 
now become straight lines for t — > oo. 



FIG. 11. (Color online) Log-log plot of the densities of active 
sites for model (A) with a = 1 and fc out = 3, for times t = 
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FIG. 10. (Color online) Log-linear plot of P SU rv(t), the proba- 
bility that there are still some active sites at time t, for model 
(A) with (j = l. All curves for k out > k c (including the critical 
one!) seem to become horizontal with P aurv (oo) > 0. 



FIG. 12. (Color online) Log-log plot of the densities of im- 
mune sites for model (A) with the same parameters as in 

Fig. [nj 



mainly depends on long range infections, as infections 
between nearest neighbors just fill in gaps in the cluster of 
immunes. It seems that r\ ~ log A: out for large fc out , so that 
n(t) ~ exp(ln/c out ln<). More precisely, both data sets 
are very well fitted by r\ = 0.9 + 0.841n(fc out - fe out , s hort), 
where fc ou t, short is the number of short-range infections. 
We have no theoretical argument for this. 

Survival probabilities P SU rv (t) for model (A) are shown 
in Fig. 10 As expected, they all tend to constants when 
the process is supercritical. More surprising, it seems 
than also for the critical case fc out = k c = 2.6 we have 
P surv (t) — ¥ const for t — > oo. Indeed, even in the clearly these distributions look similar to those for the supercrit- 



subcritical case fc ou t = 2.4, we see that P SU rv(^) first 
curves upward, before it finally goes to zero. This is 
in contrast to the behavior of SIR epidemics with short 
range contacts, but it is precisely what was proven rig- 
orously in |33| : When a = 1, the percolation transi- 
tion is discontinuous in the sense that the order parame- 
ter (which is just lim^oo P SU rv(0) jumps discontinuously 
when the control parameter fc out passes through the crit- 
ical point. Furthermore, P surv (oo) > 2/k c as predicted 
in [33] . 

Spatial distributions of active and immune sites for 
model (A) are shown in Figs. 11 and 12 Superficially, 
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ical case with a < 1 shown in Figs. [2] and [3J but there 
are important differences: 

• The curves are now equidistant for exponentially 
increasing times, corresponding to the fact that the 
number of active sites, the number of immune sites, 
and the effective radius all increase like powers of £, 
while they increased faster than polynomial when 
a<l. 

• While the peaks in Fig. [2] become narrower with 
increasing t, now the shapes of the curves are in- 
dependent of t, suggesting finite time scaling for 
(7 = 1, but not for o < 1. This is also seen by 
making formal data collapses (not shown here). 

In view of the fact that R(x, t) is constant [t- 
independent) for small values of x we make for 
a = l and fc out > fc c the ansatz 



R(x,t) * </>(x/£(t)) 
for t — > oo, with 



const for \z\ <C 1, 
Iz]-*- 1 for \z\ > 1. 



(14) 



(15) 



and £(t) — > oo for t — > oo. From this ansatz follows 
immediately that has to scale exactly like N(t), 



(16) 



and that the density of active sites satisfies the scal- 
ing 



p(i,i)«t-V(i/f(t)) 



(17) 



with ip(z) = zcj)'(z). 



• Since, according to Eqs.(17l and (|16j), both the 
peak heights and the peak positions of p(x,t) scale 
as powers of t for a = 1, also the peak heights must 
scale as a power of the positions, 



f (t) = maxp(x,t) ~ £(t)-V(l+n) 



(18) 



This is clearly seen in Fig. 11 where we obtain 
1/(1 + rj) = 0.369(3) or rj = 1.71(2), in per- 
fect agreement with the direct measurement r\ = 
1.692(6). But as we shall see later, Eq. (18) does 
not hold for a < 1. 

Since the generated clusters are basically compact, we 
can measure rj also by measuring the average of log|x|, 
either over the active or over the immune sites. Both 
averages should scale as log £(£) ~ (1 + rj) log t. This was 
indeed verified numerically. 

Equations ( 14 1 to ( 18 1 were also tested for several other 
values of fc ou t, and were satisfied in all cases. This in- 
cludes even the critical case k c . Also there, R(x, t) is flat 
for x <C £,(t), showing that the pair connectedness does 



not decrease as long as x < £(t). This is of course in 
agreement with the result that P S urv does not decrease 
to zero at the critical point, and that the percolation 
transition is discontinuous for a = 1 |33] . 



We should point out that Eqs. (14) to (18) hold also 



for model (B) if a = 1, although finite time corrections 
are somewhat larger since not all links are distributed 
according to a power law. For one typical case (fc ou t = 
3.0), see Fig. 13 In spite of the visible deviation in the 
curve for t = 20, the peaks line up for larger t along 
a perfect power law, with -q = 1 — 1/0.610(4) = 0.64(1). 
This should be compared to rj = 0.6484(4) from the direct 
measurement of n(t). 

Finally, let us discuss the predictions of mean field the- 
ory for the case a = 1. The exact evolution for R(x,t) 
can be written as 



R(x,t + 1) = R(x,t) 

+ ^2P(x - y)(n act (y,t)(l - n 



immune 



(19) 

0M))> , 



where n act (x,t) and ni mmuno (x, t) are the exact fluctu- 
ating densities of active and immune sites, and angular 
brackets indicate an ensemble average. The mean field 
assumption is 



immune 



(M)) 



(nact(y,f))(n 
p(x, t)R(x, t). 



immune 



>,*)) 



(20) 



Inserting this and the scaling ansatzes into the r.h.s. of 



Eq. ( 20 ) , we find that it gives back the scaling ansatz for 
the l.h.s., i.e. our scaling assumptions are consistent with 
mean field theory. In contrast, the scaling ansatzes would 
not be compatible with this mean field theory, if cr ^ 1. 



IV. THE CRITICAL MODEL (A) FOR < a < 1 

Critical percolation with long range infection has been 
studied both by means of field theory [9l [10] and by 
means of simulations [TD]. The field theoretic analysis 
(using the epsilon expansion) should hold in d dimension 
of space for any d, provided a sa d/3 ® [TU] - thus it 
should also be applicable to d = 1, provided a w 1/3. 
Remarkably, it predicts that mean field theory does not 
only hold for a < (where it holds for the supercritical 
case), but also for < a < 1/3. It also predicts, for 
d > 2, the value of a above which short range behav- 
ior should be observed. Again it is remarkable that this 
value is not where short range behavior is observed in 
the supercritical case (a = d), but at a smaller value of 
a [ID]. 

While these predictions have been reasonably well con- 
firmed for d = 2 |10| (much larger simulations, again on 
lattices with 2 64 sites, will be published elsewhere [35j), 
it seems that no simulations were yet done for d = 1. 

The most obvious strategy for finding the critical val- 
ues k c of fc ou t is to start again from a single seed and to 
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FIG. 13. (Color online) Log-log plot of the densities of ac- 
tive sites for model (B) with a = 1 and k out = 3, for times 
t = 20, 40, . . . 20480. Compared to Fig. [TT] the curves for the 
smallest t values now show less perfect scaling. 



FIG. 15. (Color online) Log-linear plot of t 



i(t) for 



k^t = -1 .060805 



k out = 1 .060790 




10000 



FIG. 14. (Color online) Log-linear plot of n(t) for model (A) 
with a = 1/3 and k out « k c . These data are compatible 
with k c — 1.06799(1) and r\ = 0, if we allow for logarithmic 
corrections, nit) ~ [mt] Q with a Ri 0.5 to 1. 



look for the best scaling behavior n(t) ~ t v in the large-r; 
limit. In general this works without too much problems, 
but we have to expect large finite-t corrections near any 
change of the critical behavior, i.e. in particular near 
o" = l/3 and near a = 1. But one must be aware of 
surprises. As two examples we show n(t) versus t for 
cr = l/3 (where we expect rj = 0, but slow convergence 
due to possible logarithmic corrections) and for a = 0.4, 
where we should a priori expect much less problems. The 
data, based on simulations of typically ss 10 8 clusters, are 
shown in Figs. 14 and 15 While Fig. [Ml is compatible 



model (A) with a — 0.4 and fc ou t ~ k c . These data are 
compatible with k c = 1.094538(2) and r) = 0.105(3). The 
main uncertainty of 77 comes from the obvious corrections to 
scaling. Notice, however, the expanded y-scale. Notice also 
that finite (lattice) size corrections are negligible in these sim- 
ulations, even for the largest values of t. Thus, in order to 
obtain the true critical exponents, one should take the large-t 
behavior serious, and should not try to fit in an intermediate 
region of t. 



with nothing more than the expected logarithmic correc- 
tions, the corrections to scaling seen in Fig. [15] are much 
more complicated. They cannot be described by a single 
power term, and they contribute most to the uncertainty 
of our estimate of 77. 



Plots similar to Figs. 14 and [15] were also made for 
several other values of a. In addition to n(t), we also 
looked in the same way at P surv (f ) (plots not shown) , for 
which we assumed 



-a 



(21) 



Results are shown in Figs. [16] and [17] They first of all 
confirm the prediction that mean field behavior holds for 
all cr < 1 /3, not only for a < as in the supercritical case. 
Secondly, they confirm our results of the last section, in 
particular we see that 5 — > when a — > 1. On the basis 
of Fig. [17] we conjecture more precisely that S s» 1 — a 
when a — > 1. Finally we compare to the predictions of [S], 
represented in Figs. [16] and [17] by straight lines. These 
predictions are based on an e-expansion with e = d — 3a. 
Although this expansion is most likely only asymptotic, it 
should give correct results when e — > 0. This is definitely 
the case for Fig. [16] where the prediction seems to be 
valid up to a rs 1. For (5 the agreement is much worse. 
Even for a — 0.4 it gives a prediction for 1 — S that is 
too large by « 15 per cent. This might be related to 
the anomaly seen in Fig. |15[ and might indicate that 
we have still underestimated finite-£ corrections in the 
regime 1/3 < a < 1/2. The latter is also suggested by 
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FIG. 16. (Color online) Exponents rj(a), denned via n(t) ~ t v , 
versus a for critical model (A) . The horizontal line is the mean 
field prediction r) = for a < 1/3. The straight tilted line 
is the prediction of [5] that should be exact for a slightly 
above 1/3. Finally for a — 1 we have the value 77 = 1.428(12) 
obtained in Sec. III. 



FIG. 18. (Color online) Densities of immune sites for critical 
epidemics with a — 1/2, for t — 10, 10 2 , . . . 10 4 . For large x 
all curves show the by now well known tail R oc l/x 1+cr , but 
for smaller values of x the flat part seen in Figs. 3, 4 and 11 
is replaced by a power law R oc l/x 1- " . 



1 

0.8 
0.6 
0.4 
0.2 




direct [Eq. (16)] 
exact 
e-expansion 
rhs. of Eq.(28) 



_l 1 1 1 1 l_ 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

o 



FIG. 17. (Color online) Exponent S versus a for critical model 
(A). The points with error bars are obtained directly from the 
definition P surv (t) ~ t~ s , while the dotted c urve is obtained 
from the data shown in Fig. 16 using Eq. (331. Again the 



straight lines are the exact mean field prediction (for a < 1/3) 
and the epsilon expansion result of [5] (for a slightly above 
1/3). 



the slight disagreement in the region 1/3 < a < 1/2 with 
Eq. pi 



For SIR epidemics with short range infection, one can 
also define a critical exponent z that describes how the 
correlation length increases with time, £ ~ t z . Figure 
[T8| suggests that there is indeed a well defined diverging 
length scale (the location where the break of slopes oc- 



curs; notice that this length scale should be defined by 
geometric averaging [8] , instead of the arithmetic averag- 
ing usually taken for short range infects - characterizing 
length scales by arithmetic averages when a < 1, as done 
e.g. in [Hj, can lead to dubious results). In that figure 
we present the densities of immune sites for different val- 
ues of t, for one randomly picked value of a. The data 
shown in Fig. 18 are for a — 1/2, but similar results were 
found for all other < < 1. In particular, in all cases 
we see two different powers for x < £(t) and x > £(t). 

More precisely, we observe the by now well known tail 
R(x,t) oc l/x 1+cr for large x, but for small x we observe 
a completely new phenomenon. Instead of the flat parts 
seen in Figs. 3, 4 and 11 for small x we now see another 
power law, 



R(x,t) oc 1/x 1 ' 



(22) 



for x < £(t), where is a function that increases like 
a power for large t, 



(23) 



Indeed, data collapses (not presented here) show that 

(24) 



R(x, t) satisfies a scaling law similar to Eq. 14 
R(x,t) = x^^ix/^t)) 

where 



$(ar) = ^ w ; " ^ ' (25) 
V ; \ \x\~ 2 ° for |a;| > 1. 

It is not difficult to relate the exponent z to the expo- 
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FIG. 19. (Color online) Critical exponent z versus a. As in 
Fig. |17[ points with error bars are from direct measurements 
(using now Eq. ( 28 1 ) , while the dotted curve is the prediction 
from Eq.p7l). 



ncnt rj. Indeed, by summing R{x, t) over all x we obtain 
N(t) = ^R(x,t) 



m dx , , „„„ r°° dx 



o a a- 

a a 



l+cr 



(26) 



where we used that < a < 1. Since N(t) ~ t 1+r > , we 
have thus [H] 



(27) 



FIG. 20. (Color online) Densities of active sites for critical 



epidemics with a = 1/2, for t = 10, 10 2 , 
line has the slope predicted by Eq. (30 1. 



10 4 . The straight 



from which we find that p(x, t) for fixed t has a peak at 
x « £(t) with a height 



-1-(1-ct)z 



a-l-l/z 



(30) 



This prediction is numerically verified for a = 0.5 in 
Fig.0 

Finally, we can relate the exponents t] and S by a hy- 
perscaling relation as follows: Let us consider the growth 
of two clusters, one starting at position and the other 
at x, for some large time t. If x — £(t), the probability 
that they overlap (i.e. have at least one site in common) 
is of order 1, provided both of them have survived up to 
t. Thus, up to constants of order one, 

R(x,t) ~ at) 1 '* ™ [Pit)? for x = £(t) (31) 



For a — > 1 this gives £(t) ~ t 1+v as found in Sec. Ill, 
while it gives z — > oo for a — > 0, indicating that in this 
limit £(i) is infinite and R{x, t) ~ 1/t is a pure power law 
for all t. 

As in the case a = 1 discussed in the previous section, 
we can also measure z directly by measuring the average 
logarithmic distance of active sites from the origin, 



X act (t) = exp[(p(x,t) ln\x\)]~t z 



(28) 



Values obtained in this way are shown in Fig. |19| In 
this figure we also show the values predicted by Eq. (27), 
finding perfect agreement. 



Using Eqs. (24) and (23) we obtain for the density of 
active sites 



p(x,t) 



„<7-l 



(29) 



which gives [HI ITU] 



(a-l)z 



i-2(5 



5=hl-a)z=±(l-<j)(l + V ) 
i la 



(32) 



(33) 



Like other hyperscaling relations, it only holds up to 
the critical dimension, which in the present case means 
that it holds for a > 1/3. For a < 1/3 the clusters are so 
sparse that they overlap with probability < O(l). Again, 
this prediction is satisfied for all tested values of a. 

As in other critical phenomena, there is one other inde- 
pendent critical exponent that governs the behavior for 
fcout 7^ k c but k out — > k c . Traditionally, for percolation 
with local contacts this can be cither the order parameter 
exponent (3 or the correlation length exponent v (for def- 
initions and relations between them, see e.g. [3S] or |10) . 
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FIG. 21. (Color online) Correlation time exponent u t against 
cr, obtained by data collapse for the number of infected sites. 
The straight line is the prediction if the the e-expansion [9]. 
For a < 1/3, one has the mean field prediction v t = 1. 



We use here the correlation time exponent v tl which is 
informally denned via the scaling of a characteristic time 
scale r ~ |fe out — k c \ Vt . Using this together with the FTS 
ansatz Eq. Q, we can then determine numerical values 
of v t by plotting t 5 n{t) versus (fe ou t — k c )t v for different 
trial values of y. Data collapse should occur for large t 
and finite (fc out — k c )t v , when y = l/v t . Values obtained 
in this way are shown in Fig. [2TJ We again see that the 
mean field prediction v t = 1 for a < 1/3 is satisfied, and 
that the e-expansion [5] is correct for a slightly larger 
than 1/3. For a — > 1 one finds that v t diverges, which 
is in agreement with the prediction [3(2 [311 133] than the 
transition is of BKT type for a = 1. 

As a last remark we should point out that the the- 
oretical discussion of the present section applies only to 
infinite systems. For critical phenomena with short range 
interactions, the correlation length £ describes both what 
happens in finite systems and, in infinite systems, at fi- 
nite times or finite distances from the critical point. In 
the present case, due to the definition of £ via a geometric 
average, it is less clear how it relates to the finite system 
size behavior. 



V. SUPERCRITICAL EPIDEMICS WITH 

< a < 1 

In contrast to the critical case and to the case with 
a = l, where the epidemics are described by scaling laws, 
the situation seems much less clear for supercritical epi- 
demics with < a < 1. 

In this case mean field theory predicts for both mod- 
els exponential growth of N(t) and of the spatial extent 
[2"Tll25| . The reason is very simple. For a < 1 the wave of 



infection propagates, in mean field theory, like a 'pulled' 
[2 3) front, i.e. the growth of the cluster is mainly con- 
trolled by its most advanced 'avant garde' (for a > 1, in 
contrast, we have seen that the front is 'pushed' by the 
region where the density is large). In this region ahead of 
the main front the density is very small, and thus satu- 
ration effects are negligible. Practically every attempted 
infection also succeeds, and thus the density increases 
exponentially with time as 



R 



k 1 

""out 



( for fixed large a;) 



(34) 



this 



Together with the spatial power law R ~ x 
means that an effective front position Xf mn t(t), defined 
by 



R(Xfrant (*))*) = Const , 

must advance exponentially, like |21) 



•Efront (^) 



„*/(l+<r) 



(35) 



(36) 



As we have seen in Sec. 
data. 



II, this is not supported by our 



The argument leading to Eq. ( 36 1 was criticized in [26 



who argued that the advance of the front should be linear 
in time. Unfortunately, this is not supported by the data 
either, and it is easy to see why. In |26j it was assumed 
that new infections are not successful, if they appear in 
regions with extremely small density. This would be cor- 
rect, if we had assumed cooperativity (or 'synergy') in the 
infection [U [37] . But this was not assumed in [211 ES] , 
nor was it stated explicitly in [26] . It is also not assumed 
in the present paper. 

We claim that the problem is, instead, a break-down 
of mean field theory. It is true that the front is pulled, 
and it is true that the average density, averaged over the 
entire ensemble, is small in the tail of the front. But as 
soon as a site far ahead of the previous front (i.e., in a 
region with very small density) is infected, it will generate 
its own little 'colony' and create locally a spot with high 
density. Thus, in spite of the very small average densities, 
there will always be non-negligible saturation effects due 
to the not-so-small actual densities. 

Let us define a time dependent growth rate a(t) by 



a(t) = In 



n(t+l/2) 
n(t - 1/2) 



(37) 



Growth rates for a = 0.5 and different values of k oni are 



shown in Fig. 22 According to Eq. (361 we should expect 



a(t) = a m f(t) 



In k n 



1 



(38) 



but we see that the measured a drops, for all tested values 
of fc ou t, far below this mean field value. Instead, it seems 
that a(t) decreases for large t like a power 



a(t) ~ aot 1 



(39) 
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FIG. 22. (Color online) Log-log plot of growth rates (as de- 
fined in Eq. (37l)for a = 0.5, and for different values of fc out . 
In all cases, the measure d g rowth rates are much smaller than 
those predicted by Eq. (38 1, except for very small t. 



FIG. 24. (Color online) Stretched exponential powers as a 
function of a for models (A) and (B). The continuous curves 
are the predictions of Biskup [38] (lower curve) and Emmerich 
et al. [19: (upper curve). 
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FIG. 23. (Color online) Log-log plot of re-scaled growth rates 
against re-scaled time. The straight line corresponds to 7 = 
0.49. 



where arj depends on fc out , but 7 is independent of k ou t. 
In particular for a = 0.5, we obtain 7 = 0.49(2). 
This means that for large t 



n(t) 



(40) 



with b = ao/7, e.g. for a = 0.5 we have a stretched 
exponential with exponent 0.49(2). This is particularly 
evident if we collapse these data onto a single curve 
by rescaling t and a by arbitrary functions of k auiil see 
Fig.[2| 

The value 7 = 0.49 is to be compared with Sec. II, 
where we obtained an exponent 0.59(1) by a straightfor- 



ward fit. The discrepancy between these two estimates 
results from systematic deviations from a pure stretched 
exponential, which are indeed visible in both Figs. [5] and 



22 



In both cases the curves bend downward (instead of 
being straight as for clean stretched exponentials), indi- 
cating that a naive fit overestimates the growth. Thus it 
seems likely that the present estimate obtained via the 
growth rate is more reliable. In addition, we made sim- 
ilar fits with n(t) replaced by N(t). Both should satisfy 
asymptotically the same stretched exponential, but with 
different power prefactors. The estimate of 7 via the 
growth rate seems more robust than the direct estimate 
of Sec. II. Finally, we estimated by both methods (direct 
fit & growth rate) the exponents for trial functions of the 
type f(x) ~ x a exp(6i 7 ) with various (positive ad nega- 
tive) prefactor powers. In most cases the growth rate 
method gave better results. 

For other values of c, very similar results were found. 
In particular, the exponent 7 is also independent of fc ou t 
within the estimated errors [35] , suggesting that it really 
depends only on a. The values obtained in this way for 
both models are shown in Fig. 24 We see that both 
models give very similar results. Both are consistent, 
within the errors, with 7 = 1 — cr, but it seems that 
this is not the correct behavior. A percolation model 
very similar to the static part of model (A) was studied 
rigorously by Biskup [38l [40] . He obtained results for the 
graph diameter of clusters embedded in large but finite 
lattices [30] and for the average graph distance between 
nodes with large Euclidean distance on infinite lattices 
[38] , Neither of this is precisely what is measured in the 
present work, but his results strongly suggest that 
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(41) 



1+<T 
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This prediction is also shown in Fig. 24 It agrees rea- 
sonably well with our simulations. The simulations are 
systematically too high, indicating that we have still sub- 
stantial finite cluster size corrections. Notice that naive 
fits like Fig. [3] would give even larger estimates. 

The graph diameters of graphs embedded in finite lat- 
tices were also measured by pll] [4XJ . Unfortunately, these 
authors used a rather complicated algorithm which pre- 
vented them from using very large lattices and from ob- 
taining high statistics. It also introduced particularly 
large finite lattice corrections, and the data were ana- 
lyzed by fitting stretched exponentials directly via plots 
like Fig. [5] It is presumably for these reasons that [T51I4T] 
obtain 7 = 1 — a 2 (in our notation), which is clearly 
incompatible with our data and with the prediction of 



ESjUI (see Fig. [24]) 02]. 

Figure [3] suggests a scaling ansatz 

R{x,t) = <j>{x/i{t)) 



(42) 



similar to Eq. ( 14 ) , but with £ (t) being a stretched expo- 
nential, 



^(t) ~ N(t) ~ exp(fcT) , 



(43) 



instead of a power law. From this ansatz follows 

p(x,t)=F- 1 il>(x/t(t)) ( 44 ) 

where ip(z) — z <j>'{z). This means in particular that 
Pmax(t), defined in Eq. (18), decreases like a power of t, 

t (t) = f- 1 . 



Pn 



(45) 



This is easy to understand. As we pointed out already 
at the beginning of this section, the effective growth rate 
decays as f 1 , because new attempted infections are in- 
creasingly more likely to target sites that are no longer 
susceptible. But there it applied to the 'front of the 
front', while here we see that it applies also to the 'core 
of the front' where most of the mass growth occurs. We 
can interpret Eq. ( 45 ) also as a manifestation of a weak 
sort of 'self-organized criticality' [33] in the sense that 
the speed of growth is such that the density at the front 
x « £(i) converges exactly to its critical value. 

Equations ( 45 1 and ( 43 ) together imply the remarkable 
relation 



Pmax(^) 



(46) 



which should hold for any value of er € [0, 1) and does 
not involve a or 7 explicitly. 

Equations ( 42 ) to ( 46 ) are tested numerically, again for 
a = 0.5 and fe out = 2, in Figs. [25] to [27] In Fig. [25] we 
see a perfect data collapse for sufficiently large t, while 



we see in Fig. 26 that 



r 



(47) 



with /3 = 0.037 seems to fit the large time asymptotics. 
Together with Eq. (43), this would however imply that 
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FIG. 25. (Color online) Collapse plot for the data (for t > 15 
only) shown in Fig. [3] We see a perfect data collapse for 
t > 50. 




FIG. 26. (Color online) Enlargement of part of Fig. [2] It 
shows that the peak heights decrease with t, if t > 20. For 
very large t, where Fig. [25] shows a data collapse, this de- 
crease is compatible with a power law p max ~ with very 
small /3, but we shall argue that this dependence is actually 
logarithmic. 



Pmax(^) decreases with t faster than a power, which is in 
compatible with Eq. (45). Therefore we show in Fig. 27 
the ratio 



,> as a function of t, in order to test 
Eq. (46). Here, £(t) is as obtained in the collapse shown 
in Fig. 25 This defines £(t) only up to a constant. This 



constant (used, by the way, also in Fig. 25) is fixed so 
that a straight line is obtained in Fig. |27| We see that 
the data follow indeed a nice linear relationship, showing 
that at least the entire scheme is internally consistent. 
Finally, before leaving this section, let us point out that 
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FIG. 27. (Color online) Log-log plot of p ma x,ln£, and of 
ln£/r/io max against t. According to Eq. (46 1 the latter should 



be proportional to t in the large t limit, which according to 
Fig. [25] should be reached for t ~ 50. The straight line has 
slope 1. 
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FIG. 28. (Color online) Log-log plot of R(X, i) versus x, 
for t = 10,20,50, 100, 200, 350, 550, 800, 1000, 1200,and 1400. 
Control parameters are a — 0.25 and k out = 0.04. The latter 
is ca 20% above fe c . 



the scaling relation (|42j) sets in very late when k out w k c , 
in particular when a < 1/3. In this range of a, mean field 
theory hold for fc out = k c , but it does not hold for any 
&out > k c . Thus, the cross-over from mean-field to 'semi- 
local' behavior has to happen in a very narrow region 
with fc out — k c <C 1. We illustrate this with Fig. [28} 
where we show R{x, t) for a = 0.25 and fc ou t = 0.04 
versus x for different values of t, in a way completely 
analogous to Fig. 3 (where a = 0.5, A; out = 2). Instead 
of the structureless curves in Fig. 3, we see now that for 



short times the system behaves as if it were critical and 
mean field. For t > 200 it becomes clear that the process 
is supercritical, but newly infected sites are very far from 
the seed. It is only for t > 1000 that the density starts 
to grow again appreciably at intermediate distances x ~ 
10 6 , in order to reach finally its asymptotic value O(l). 



VI. DISCUSSION AND CONCLUSIONS 

Networks embedded in space with connections which 
preferentially link close neighbors but have also non- 
vanishing chances to link nodes far apart have numerous 
applications, from biology to social sciences. Very early 
it was already proposed to model the link length distribu- 
tion by power laws [8| I44j . In the present paper we gener- 
ate such networks by epidemic processes with power be- 
haved distributions for contacts, i.e. infections occur over 
distances x whose probability decays as P(x) ~ l/x 1+a . 

While the case of two spatial dimensions will be treated 
in a forthcoming paper [35] , we restricted ourselves here 
to d = 1. This is of course less interesting from the point 
of view of applications, but it allows much more detailed 
and precise analysis. Our results concern mainly three 
different regimes: 

(i) critical epidemics with < a < 1, 

(ii) supercritical epidemics with a = 1, 

(iii) supercritical epidemics with < a < 1. 

In all three cases we obtain significant new results, by 
using simulations on unprecedentedly large sizes. Using 
hashing, the lattices we used have 2 64 « 1.8 x 10 19 sites. 
Such large lattices are needed in order to avoid finite size 
effects, in view of infections that can spread over billions 
of sites in one time step. 

For critical epidemics we verify predictions from field 
theory and show that there is one independent criti- 
cal exponent less than for critical epidemics with short 
range contacts in > 2 dimensions (in d = 1, epidemics 
with short range contacts die out). For supercritical 
epidemics with a = 1 we verify old predictions based 
on the Fortuin-Kasteleyn connection between percolation 
and Potts models. And for supercritical epidemics with 
< a < 1 we verify predictions [381 HO] that they lead 
to 'medium-size world' networks, i.e. to networks that 
are not 'small world' in the sense that their size grows 
exponentially with their graph diameter, but which are 
also not fractal in the sense that this mass grows like 
a power. Instead, it grows like a stretched exponential. 
Related to this is the observation that supercritical epi- 
demics with < a < 1 spread neither with fixed nor 
with exponentially increasing velocity, in contradiction 
to previous claims. 

In all three cases we find several new scaling laws that 
are strongly suggested numerically and which we show 
in some cases to satisfy non-trivial consistency relations, 
but for which we do not give theoretical derivations. The 
need to provide these proofs is one of the main open 
problems. 



15 



Finally, a last problem which we left open is finite size 
behavior. Our strategy was to use lattice sizes which 
are big enough so that finite (lattice-)size effects can be 
safely neglected. In view of the interplay between the 
length scales set by the lattice size, the critical correla- 
tion length, and the large contact distances, we expect 
finite size effects to be not as simple as in conventional 
finite size scaling for critical phenomena with short range 
interactions. 
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